Charged polytropic compact stars. 
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In this work, we analyze the effect of charge in compact stars considering the limit of the maximum amount 
of charge they can hold. We find that the global balance of the forces allows a huge charge (~ 10 20 Coulomb) 
to be present in a neutron star producing a very high electric field (~ 10 21 V/m). We have studied the particular 
case of a polytropic equation of state and assumed that the charge distribution is proportional to the mass density. 
The charged stars have large mass and radius as we should expect due to the effect of the repulsive Coulomb 
force with the M/R ratio increasing with charge. In the limit of the maximum charge the mass goes up to ~ 10 
Mq which is much higher than the maximum mass allowed for a neutral compact star. However, the local effect 
of the forces experienced by a single charged particle, makes it to discharge quickly. This creates a global force 
imbalance and the system collapses to a charged black hole. 
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I. INTRODUCTION 

In 1924, Rosseland| 1] studied the possibility of a self gravi- 
tating star on Eddington's theory to contain a net charge where 
the star is modeled by a ball of hot ionized gas (see also Ed- 
dington |2]). In such a system the electrons (lighter particles) 
tend to rise to the top because of the difference in the par- 
tial pressure of electrons compared to that of ions (heavier 
particles). The motion of electrons to the top and further es- 
cape from the star is stopped by the electric field created by 
the charge separation. The equilibrium is attained after some 
amount of electrons escaped leaving behind an electrified star 
whose net positive charge is of about 100 Coulomb per so- 
lar mass, and building an interstellar gas with a net negative 
charge. As shown by Bally and Harrison j3|, this result ap- 
plies to any bound system whose size is smaller than the De- 
bye length of the surrounding media. The conclusion is that 
a star formed by an initially neutral gas cannot acquire a net 
electric charge larger than about 100C per solar mass. It is 
expected that the sun holds some amount of net charge due 
to the much more frequent escape of electrons than that of 
protons. Moreover, it is also expected that the escape would 
stop when the electrostatic energy of an electron e<I> is of the 
order of its thermal energy kT. This gives for a ball of hot 
matter with the sun radius, a net charge Q ~ 6.7 x 10 _6 T 
(in Coulomb). Hence, the escape effect cannot lead to a net 
electric charge much larger than a few hundred Coulomb for 
most of the gaseous stars. 
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For Newtonian stars, the net charge of 100 C per solar mass 
is obtained by the balance between the electrostatic energy 
eQ/r and the gravitational energy mM/r (Glendenning|4|). 
However, for very compact stars, the high density and the rel- 
ativistic effects must be taken into account |5]. In a strong 
gravitational field, the general relativistic effects are felt and 
the star needs more charge to be in equilibrium. Moreover, for 
very compact stars, the induced electric field can be substan- 
tially higher than in the case of the sun. For instance, the same 
amount of charge yields an electric field approximately 10 9 
times larger at the surface of a neutron star than at the surface 
of the sun. So, even a relatively small amount of net charge on 
compact stars can induce intense electric fields whose effects 
may become important to the structure of the star. This fact 
deserves further investigation. 

The general relativistic analog for charged dust stars was 
discovered by Majumdar |6] and by Papapetrou 0|, and fur- 
ther discussed by Bonnor |8] and several other authors|9]. 
Study for the stability of charged fluid spheres have been done 
by BekensteinU, Zha ng e t al.fioll. de Felice & Willi. Yu & 
LiuO, de Felice et al.fOl, Anninos & Rothm anll4ll and oth- 
ers. This was indirectly verified by Zhang et al. 1 10] who found 
that the structure of a neutron star, for a degenerate relativistic 
fermi gas, is significantly affected by the electric charge just 
when the charge density is close to the mass density (in ge- 
ometric units). In the investigations by de Felice et al., and 
by Anninos & Rothman, they assumed that the charge distri- 
bution followed particular functions of the radial coordinate, 
and they were mostly interested in the extreme Q = M case. 

Our basic consideration to incorporate charge into the sys- 
tem is in the form of trapped charged particles where the 
charge goes with the positive value. The effect of charge does 
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not depend on its sign by our formulation. The energy den- 
sity which appears from the electrostatic field will add up to 
the total energy density of the system, which in turn will help 
in the gaining of the total mass of the system. The modi- 
fied Tolman-Oppenheimer-Volkoff (TOV) equation now has 
extra terms due to the presence of the Maxwell-Einstein stress 
tensor. We solve the modified TOV equation for polytropic 
equation of state (EOS) assuming that the charge density goes 
with the matter density and discuss the results. The formation 
of this extra charge inside the star is however left open. A 
mechanism to generate charge asymmetry for charged black 
holes has been suggested recently by Mosquera Cuesta et al. 
flUl and the same may be applied for compact stars too . 

This article is arranged in the following way. In SectionllTI 
we show the basic formalism for the modified TOV. In Section 
[TTTl we used this modified TOV on a polytropic EOS, discuss 
the results and the stability of the charged stars. Finally we 
make our conclusions in Section llVl 



II. THE MODIFIED HYDROSTATIC EQUILIBRIUM 
EQUATION 

We take the metric for our static spherical star as 



is the electric field. So, the total charge of the system is 
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where R is the radius of the star. 

The mass of the star is now due to the total contribution of 
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the energy density of the matter and the electric energy (^) 
density. The mass takes the new form as 



Mtot(r) = I 4*7* ( - 
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and the metric coefficient is given by 
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The stress tensor is conserved (T^ = 0). Hence, one gets 
the form of the hydrostatic equation from it as : 



dP 
dr 
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ds 2 = e u c 2 dt 2 - e x dr 2 - r 2 {d9 2 
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The stress tensor T{? will include the terms from the 
Maxwell's equation and the complete form of the Einstein- 
Maxwell stress tensor will be : 

TS = (P + eKu. + P^ + i (V lQ F Q „ - -S^F^F^ 

(2) 

where P is the pressure, e is the energy density (=pc 2 ) and u-s 
are the 4-velocity vectors. For the time component, one easily 
sees that u t — e~ v l 2 and hence u*it t = —1. Consequently, 
the other components (radial and spherical) of the four vector 
are absent. 

Now, the electromagnetic field is taken from the Maxwell's 
field equations and hence they will follow the relation 
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where is the four-current density. Since the present choice 
of the electromagnetic field is only due to charge, we have 
only F 01 = — F 10 , and the other terms are absent. In gen- 
eral, we can derive the electromagnetic field tensor F^ v from 
the four-potential A^. So, for non vanishing field tensor, the 
surviving potential is Aq = <fi. We also considered that the 
potential has a spherical symmetry, i.e., tj> = (f>(r). 

The nonvanishing term in Eq.(|3} is when v=x. This gives 
the electric field for both the t and r components as : 
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where, 



U{r) = \ [ ^r 2 p ch e x l 2 dr. 
r l Jo 



(4) 



We solve the Eqns. ||6][7]&[8J simultaneously to get our results 
for the charged compact stars. 



III. EFFECT OF CHARGE ON POLYTROPIC STARS 

We study here the effect of charge on a model independent 
polytropic EOS. We assume the charge is proportional to the 
mass density (e) like 



Pck = /xe 



(9) 



where e = pc 2 is in [MeV/fm 3 ]. In geometrical units, this can 
be written as 



p ch = ax p 



(10) 



where charge is expressed in units of mass and charge den- 
sity in units of mass density. This a is related to our charge 
fraction / as 

a = fx °- 22 ^ 36 =/x 0.86924 xlO 3 . (11) 



G 

Our choice of charge distribution is a reasonable assumption 
in the sense that large mass can hold large amount of charge. 
The polytropic EOS is given by 



P = K P 1 + 1 ' 



(12) 



where n is the polytropic index and is related to the exponent 
r as T = 1 + — . In the relativistic regime, the allowed value of 
r is | to | . We have considered the adiabatic case of T = | 
and the corresponding value of n is 1.5. Primarily, our units 
of matter density and pressure are in MeV/fm 3 . We chose a 
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FIG. 1: Central density against mass for different values of the factor 
/• 



value of k as 0.05 [fm] 8 / 3 for our polytrope that reproduces 
quite well realistic EOS for neutron stars 11611 . It should be 
noted that the amount of charge we find implies a very small 
ratio Z/A ~ 10~ 18 which justifies to use an EOS which is 
calculated for neutral matter. 

In Fig.Q, we plot the mass as function of the central den- 
sity, for different values of the charge fraction /. For the 
charge fraction / = 0.0001, we do not see any departure 
on the stellar structure from that of the chargeless case. This 
value of / is critical because any increase in the value beyond 
this, shows enormous effect on the structure. The increase of 
the maximum mass of the star is very much non-linear, as can 
be seen from the Fig. Q. 




FIG. 2: Mass as a function of radius, for different values of the factor 
/• 

In Fig.© we plotted the mass-radius relation. Due to the 
effect of the repulsive force, the charged stars have large ra- 
dius and larger mass as we should expect. Even if the ra- 
dius is increasing with the mass, the M/R ratio is also increas- 
ing, but much slower. For the lower charge fractions, this in- 
crease in the radius is very small, but a look at the structure 
for the fraction / = 0.001 reveals that for a mass of 4.3 M©, 
the radius goes as high as 35 km. Though the compactness 
of the stars are retained, they are now better to be called as 
charged compact stars rather than charged neutron stars. 
The charge fraction in the limiting case of maximally allowed 
value goes up to / = 0.0011, for which the maximum mass 
stable star forms at a lower central density even smaller than 
the nuclear matter density. This extreme case is not shown in 



Fig. (0 because the radius of the star and its mass is very high 
(68 km & 9.7 M Q respectively) which suppresses the curves 
of the lower charge fractions due to scaling. For this star, the 
mass contribution from the electric energy density is 10% than 
that from the mass density. It can be checked by using rela- 
tion (lilt that this charge fraction / = 0.0011 corresponds to 
p c h = 0.95616 x p in geometrical units. 
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FIG. 3: The variation of the charge with mass for different /. 

The Qx M diagram of Fig.Q shows the mass of the stars 
against their surface charge. We have made the charge den- 
sity proportional to the energy density and so it was expected 
that the charge, which is a volume integral of the charge den- 
sity, will go in the same way as the mass, which is also 
a volume integral over the mass density. The slope of the 
curves comes from the different charge fractions. The nature 
of the curves in fact reflects that charge varies with mass 
(with the turning back of the curves all falling in the 'un- 
stable zone' and is not taken into consideration). If we con- 
sider that the maximum allowed charge estimated by the con- 
dition (U ~ VSnP < \/87re) for ^ to be negative (Eq. ©), 
we see that the curve for the maximum charge in Fig.@ 
has a slope of 1:1 (in a charge scale of 10 20 Coulombs 1 ). 
This scale can easily be understood if we write the charge 
as Q = VGM q1 ^L ~ lO 20 ^ Coulombs. This charge Q 
is the charge at the surface of the star where the pressure and 
also are zero. So, at the surface, the Coulomb force is es- 

ar ' 

sentially balanced by the gravitational force and the relation 
of the charge and mass distribution we found is exactly the 
same for the case of charged dust sphere discussed earlier by 
Papapetrou|7| and Bonnor|8|. 

The total mass of the system M to t increases with increasing 
charge because the electric energy density adds on to the mass 
energy density. This change in the mass is low for smaller 
charge fraction and going up to 7 times the value of chargeless 
case for maximum allowed charge fraction / = 0.0011. The 
most effective term in Eq.ljSJ is the factor (M to t + 4.7rr 3 P*). 
P* = P — ^ is the effective pressure of the system because 



As electric energy density and pressure needs to be of the same order, 

2 , 

so from the fine structure constant a = ^ = j^f, we get a rela- 
tion for the charge and MeV/fm 3 which comes out as 1C ~ 0.75 X 

W^lMeVfm} 1 / 2 . 
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the effect of charge decreases the outward fluid pressure, neg- 
ative in sign to the inward gravitational pressure. With the 
increase of charge, the value of P* decreases, and hence the 
gravitational negative part of Eq.(|8} decreases. So, with the 
softening of the pressure gradient, the system allows more ra- 
dius for the star until it reaches the surface where the pressure 
(and 4?) goes to zero. We should stress that because ^ can- 
not be too much larger than the pressure in order to maintain 
4^ negative as discussed before, so we have a limit on the 
charge, which comes from the relativistic effects of the gravi- 
tational force and not just only from the repulsive Coulombian 
part. 
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FIG. 4: The positive Coulomb part and the gravitational negative part 
of the pressure gradient together with the total (4fr) are shown here 
for two different values of the charge factor /. For / = 0.0005 
and 0.0008 coming from the matter part are denoted as dPg5 and 
dPg8 respectively, those from Coulomb part are dPc5 and dPc8 
respectively. The corresponding totals are dP5 and dP8. 

This effect is shown in Fig.0 where we have plotted both 
the positive Coulomb part and the negative matter part of the 
pressure gradient. The plots are for two values of the charge 
fraction / = 0.0005 and / = 0.0008. The positive part of ^ 
maintains its almost constant value because the charge frac- 
tion / is the controller of the same, and in our case, they differ 
by a very small percentage. In the negative part, the changes 
are drastic and are mainly brought by the effective pressure as 
we already discussed. 

In our high density system, the gravitational and the 
Coulomb forces are highly coupled. Although it is difficult 
to disentangle the forces, but to a common belief, it can be 



considered that the charged particles, due to their self created 
huge field, will leave the star very soon. This process will 
however lead to an imbalance of the global forces acting on 
the star, which were previously balanced by the Gravitational 
and the Coulomb forces. This process will help in the star to 
further collapse to a charged black hole. We say it charged be- 
cause, by the time the system collapses, all the charge has not 
left the star, and they get trapped inside the black hole 11711 . 



IV. CONCLUSIONS 

In our study, we have shown that a high density system 
like a neutron star can hold huge charge of the order of 10 20 
Coulomb considering the global balance of forces. With the 
increase of charge, the maximum mass of the star recedes back 
to a lower density regime. The stellar mass also increases 
rapidly in the critical limit of the maximum charge content, 
the systems can hold. The radius also increases accordingly, 
however keeping the M/R ratio increasing with charge. The 
increase in mass is primarily brought in by the softening of the 
pressure gradient due to the presence of a Coulombian term 
coupled with the Gravitational matter part. Another intrinsic 
increase in the mass term comes through the addition of the 
electric energy density to the mass density of the system. 

The inside electric field of the charged stars are very high 
and crosses the critical field limit for pair creation (Beken- 
stein |5]). However, this issue is debatable because the critical 
field has been calculated for vacuum and one does not really 
know what the value will be in a high density system. The 
stability of the charged stars are however ruled out from the 
consideration of forces acting on individual charged particles. 
They face enormous radial repulsive force and leave the star 
in a very short time. This creates an imbalance of forces and 
the gravitational force overwhelms the repulsive Coulomb and 
fluid pressure forces and the star collapses to a charged black 
hole. 

Finally, these charged stars are supposed to be very short 
lived, and are the intermediate state between a supernova col- 
lapse and charged black holes. 
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